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1.  Statement  of  the  Problem  Studied 

Functional  differential  equations  (FDE)  with  delay  provide  a  math¬ 
ematical  model  for  a  physical  or  biological  system  in  which  the  rate 
of  change  of  the  system  depends  upon  its  past  history.  Although  the 
general  theory  and  basic  results  for  FDE  have  by  now  been  thoroughly 
investigated,  the  literature  devoted  to  this  area  of  research  continues 
to  grow  very  rapidly.  The  number  of  interesting  works  is  very  large, 
so  that  our  knowledge  of  FDE  has  been  substantially  enlarged  in  re¬ 
cent  years.  Naturally,  new  important  problems  and  directions  arise 
continually  in  this  intensively  developing  field. 

The  report  summarizes  the  results  in  the  study  and  addresses  the 
need  for  further  investigation  of  generalized  solutions  to  broad  classes 
of  FDE.  The  project  concentrated  on  differential  equations  with  piece- 
wise  continuous  arguments  (EPCA),  the  exploration  of  which  has  been 
initiated  in  our  papers  a  few  years  ago.  These  equations  arise  in  an 
attempt  to  extend  the  theory  of  FDE  with  continuous  arguments  to  dif¬ 
ferential  equations  with  discontinuous  arguments.  This  task  is  also  of 
considerable  applied  interest  since  EPCA  include,  as  particular  cases, 
impulsive  and  loaded  equations  of  control  theory  and  are  similar  to 
those  found  in  some  biomedical  models.  A  typical  EPCA  contains  ar¬ 
guments  that  are  constant  on  certain  intervals.  A  solution  is  defined 
as  a  continuous,  sectionally  smooth  function  that  satisfies  the  equation 
within  these  intervals.  Continuity  of  a  solution  at  a  point  joining  any 
two  consecutive  intervals  leads  to  recursion  relations  for  the  solution 
at  such  points.  Hence,  the  solutions  are  determined  by  a  finite  set  of 
initial  data,  rather  than  by  an  initial  function  as  in  the  case  of  general 
FDE.  Therefore,  underlying  each  EPCA  is  a  dynamical  system  gov¬ 
erned  by  a  difference  equation  of  a  discrete  argument  which  describes 
its  stability,  oscillation,  and  periodic  properties.  It  is  not  surprising 
then  that  recent  work  on  EPCA  has  caused  a  new  surge  in  the  study  of 
difference  equations.  Of  significant  interest  is  the  exploration  of  partial 
differential  equations  (PDE)  with  piecewise  continuous  delays.  Bound¬ 
ary  and  initial-value  problems  for  some  EPCA  with  partial  derivatives 
were  considered  and  the  behavior  of  their  solutions  investigated.  The 
results  were  also  extended  to  equations  with  positive  definite  operators 
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in  Hilbert  spaces.  This  topic  is  of  great  theoretical,  computational, 
and  applied  value  since  it  opens  the  possibility  of  approximating  com¬ 
plicated  problems  of  mathematical  physics  by  simpler  EPCA. 

It  is  well  known  that  profound  and  close  links  exist  between  func¬ 
tional  and  functional  differential  equations.  Thus  the  study  of  the  first 
often  enables  one  to  predict  properties  of  differential  equations  of  neu¬ 
tral  type.  On  the  other  hand,  some  methods  for  the  latter  in  the  special 
case  when  the  argument  deviation  vanishes  at  individual  points  have 
been  used  to  investigate  functional  equations.  Functional  equations 
are  directly  related  to  difference  equations  of  a  discrete  argument,  and 
bordering  on  difference  equations  are  impulsive  FDE  with  impacts  and 
switching  and  loaded  equations  (that  is,  those  including  values  of  the 
unknown  solution  for  given  constant  values  of  the  argument).  The  ar¬ 
gument  deviations  of  the  EPCA  considered  in  the  project  vanish  at 
countable  sets  of  points,  and  it  would  be  interesting  to  investigate  the 
relationship  between  EPCA  and  functional  equations.  Another  deserv¬ 
ing  direction  of  future  research  is  the  exploration  of  hybrid  systems 
consisting  of  EPCA  and  functional  equations.  Futhermore,  EPCA  are 
intrinsically  closer  to  difference  rather  than  to  differential  equations. 
Equations  with  piecewise  constant  delay  can  be  used  to  approximate 
differential  equations  that  contain  discrete  delays.  It  would  be  useful  to 
draw  a  detailed  comparison  of  the  qualitative  and  asymptotic  properties 
of  differential  equations  with  continuous  arguments  and  their  EPCA 
approximations,  which  has  been  widely  used  for  ordinary  differential 
equations  and  their  difference  approximations.  Since  the  arguments  of 
an  EPCA  have  intervals  of  constancy  we  must  relinquish  smoothness 
of  the  solutions,  but  we  still  retain  their  continuity.  This  enables  us  to 
derive  a  homogeneous  difference  equation  for  the  values  of  a  solution 
at  the  endpoints  of  the  intervals  of  constancy  and  to  employ  it  in  the 
study  of  the  original  EPCA,  thus  revealing  remarkable  asymptotic,  os¬ 
cillatory,  and  periodic  properties  of  this  type  of  FDE.  Of  course,  it  is 
possible  to  further  generalize  the  definition  of  a  solution  for  an  EPCA, 
by  abandoning  its  continuity,  and  to  include  in  the  framework  of  EPCA 
the  impulsive  functional  differential  equations. 

Along  with  mildly  weakened  solutions  of  EPCA,  the  project  also  ex- 
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plored  generalized-function  solutions  of  ordinary  differential  and  func¬ 
tional  differential  equations.  The  unifying  theme  of  the  project  is  the 
development  of  theoretically  meaningful  and  potentially  applicable  gen¬ 
eralized  concepts  of  solutions  for  important  classes  of  FDE.  A  common 
feature  of  these  equations  is  that  their  arguments  have  a  fixed  point. 
Thus,  the  argument  of  a  typical  EPCA  is  the  greatest-integer  func¬ 
tion,  and  we  also  focus  on  FDE  with  linearly  transformed  arguments. 
Hence,  it  is  natural  to  pose  the  initial-value  problem  for  such  equations 
not  on  an  interval  but  at  a  number  of  individual  points.  Contrary  to 
general  functional  differential  equations,  EPCA  of  all  types  (retarded, 
advanced,  mixed,  neutral)  have  two-sided  solutions,  and  FDE  with  lin¬ 
early  transformed  arguments  possess,  under  certain  conditions,  analytic 
or  entire  solutions.  Some  methods  in  the  theory  of  entire  solutions  are 
applied  to  prove  stability  theorems  for  linear  EPCA  with  variable  co¬ 
efficients.  Integral  transformations  establish  close  connections  between 
entire  and  generalized  functions  (distributions).  Therefore,  a  unified 
approach  may  be  used  in  the  study  of  both  distributional  and  entire 
solutions  to  some  classes  of  linear  ordinary  and  functional  differential 
equations. 

Recently  there  has  been  considerable  interest  in  problems  concern¬ 
ing  the  existence  of  solutions  to  differential  and  functional  differential 
equations  in  various  spaces  of  generalized  functions.  Many  important 
areas  in  mathematics  and  theoretical  physics  employ  the  methods  of 
distribution  theory.  Generalized  functions  are  continuous  linear  func¬ 
tionals  on  spaces  of  infinitely  smooth  functions  with  certain  conditions 
of  decay  at  infinity.  They  provide  a  suitable  framework  where  major 
analytical  operations  such  as  differentiation  can  be  performed.  Fur¬ 
thermore,  the  importance  of  the  class  of  generalized  functions  stems 
from  the  fact  that  it  includes  the  set  of  regular  distributions  repre¬ 
sented  by  locally  integrable  functions.  There  is  an  abundance  of  sin¬ 
gular  distributions,  and  the  Dirac  delta  function  is  one  of  them.  It  is 
well  known  that  normal  linear  homogeneous  systems  of  ordinary  dif¬ 
ferential  equations  (ODE)  with  infinitely  smooth  coefficients  have  no 
singular  distributional  solutions.  However,  these  solutions  may  appear 
in  the  case  of  equations  whose  coefficients  have  singularities.  We  de- 
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velop  the  methods  of  study  and  establish  some  major  results  for  linear 
ODE  in  the  space  of  finite-order  distributions  (finite  linear  combina¬ 
tions  of  the  delta  function  and  its  derivatives).  An  existence  criterion 
of  such  solutions  for  any  linear  ODE  is  found.  Necessary  and  sufficient 
conditions  are  discovered  for  the  simultaneous  existence  of  solutions  to 
linear  ODE  in  the  form  of  rational  functions  and  finite-order  distribu¬ 
tions.  The  results  are  also  used  in  the  study  of  polynomial  solutions 
to  some  important  classical  equations.  Then  distributional  solutions 
of  certain  classes  of  ODE  and  FDE  are  presented  as  infinite  series  of 
the  delta  function  and  its  derivatives.  Existence  and  nonexistence  the¬ 
orems  in  spaces  of  infinite-order  distributions  are  obtained  for  linear 
equations  with  polynomial  coefficients  and  used  to  explore  their  entire 
solutions.  We  emphasize  and  investigate  the  conditions  when  linear 
FDE  with  polynomial  coefficients  and  linearly  transformed  arguments 
have  entire  solutions  of  zero  order.  This  is  a  remarkable  dissimilar¬ 
ity  between  the  behavior  of  FDE  and  ODE  since  first-order  algebraic 
ODE  have  no  entire  transcendental  solutions  of  order  less  than  | .  An 
equally  striking  phenomenon  is  the  existence  of  distributional  solutions 
for  linear  homogeneous  FDE  without  singularities  in  the  coefficients.  In 
other  words,  distributional  solutions  to  linear  homogeneous  FDE  may 
be  originated  either  by  singularities  of  their  coefficients  or  by  argument 
deviations.  Recent  studies  have  shown  that  nonexistence  of  infinite  or¬ 
der  distributional  solutions  for  linear  time-dependent  delay  equations 
with  real  analytic  coefficients  implies  nonexistence  of  small  solutions 
(approaching  zero  faster  than  any  exponential  as  t  tends  to  infinity), 
which  is  important  in  the  qualitative  theory  of  FDE.  It  would  be  nice 
to  extend  the  results  on  distributional  and  entire  solutions  and  their 
interplay  to  partial  differential  equations  (PDE).  As  a  first  step  in  this 
direction,  one  could  take  a  linear  PDE  in  two  independent  variables 
with  polynomial  coefficients  that  admits  separation  of  variables,  then 
consider  a  series  whose  terms  are  products  of  distributional  solutions 
of  the  ordinary  differential  equations  arising  after  separation. 
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2.  Summary  of  the  Most  Important  Results 

We  shall  describe  now  some  of  the  work  that  has  been  done  in  the 
project  on  the  differential  equations  that  we  call  equations  with  piece- 
wise  continuous  arguments,  or  EPCA.  A  brief  survey  of  the  present 
status  of  this  research  is  given  in  [5] 

A  typical  EPCA  is  of  the  form 

x'{t)  =  f(t,x(t),x(h(t))), 

where  the  argument  h(t)  has  intervals  of  constancy.  For  example,  in  [4] 
equations  with  h(t )  =  [<],  [<  —  n],  t  —  n[t]  were  investigated,  where  n 
is  a  positive  integer  and  [•]  denotes  the  greatest-integer  function.  Note 
that  h(t)  is  discontinuous  in  these  cases,  and  although  the  equation  fits 
within  the  general  paradigm  of  delay  differential  or  functional  differen¬ 
tial  equations,  the  delays  are  discontinuous  functions.  Also  note  that 
the  equation  is  nonautonomous,  since  the  delays  vary  with  t.  Moreover, 
as  we  have  mentioned,  the  solutions  are  determined  by  a  finite  set  of 
initial  data,  rather  than  by  an  initial  function,  as  in  the  case  of  gen¬ 
eral  FDE.  In  fact,  EPCA  have  the  structure  of  continuous  dynamical 
systems  within  intervals  of  certain  lengths.  Continuity  of  a  solution 
at  a  point  joining  any  two  consecutive  intervals  then  implies  recursion 
relations  for  the  solution  at  such  points.  Therefore,  EPCA  represent  a 
hybrid  of  continuous  and  discrete  dynamical  systems  and  combine  the 
properties  of  both  differential  and  difference  equations. 

An  equation  in  which  x'(t )  is  given  by  a  function  of  x  evaluated  at  t 
and  at  arguments  [<],...,  [t  —  IV],  where  N  is  a  non-negative  integer, 
may  be  called  of  retarded  or  delay  type.  If  the  arguments  are  t  and 
[i-f  1], . . . ,  [t+ N],  the  equation  is  of  advanced  type.  If  both  these  types 
of  arguments  appear  in  the  equation,  it  may  be  called  of  mixed  type. 
If  the  derivative  of  highest  order  appears  at  t  and  at  another  point, 
the  equation  is  generally  said  to  be  of  neutral  type.  All  types  of  EPCA 
share  similar  characteristics.  First  of  all,  it  is  natural  to  pose  the  initial- 
value  problem  for  such  equations  not  on  an  interval  but  at  a  number 
of  individual  points.  Secondly,  for  ordinary  differential  equations  with 
a  continuous  vector  field  the  solution  exists  to  the  right  and  left  of  the 
initial  t- value.  For  retarded  FDE,  this  is  not  necessarily  the  case  [6]. 
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Furthermore,  it  appears  that  advanced  equations,  in  general,  lose  their 
margin  of  smoothness,  and  the  method  of  successive  integration  shows 
that  after  several  steps  to  the  right  from  the  initial  interval  the  solution 
may  even  not  exist.  However,  two-sided  solutions  do  exist  for  all  types 
of  EPCA.  Finally,  the  problems  for  EPCA  studied  so  far  are  closely 
related  to  ordinary  difference  equations  and  indeed  have  stimulated 
new  work  on  these. 

It  is  important  to  note  that  EPCA  provide  the  simplest  examples 
of  differential  equations  capable  of  displaying  chaotic  behavior.  For 
instance,  following  Ladas  [7],  one  can  see  that  the  unique  solution  of 
the  initial-value  problem 

x'(t)  =  3z([*])  -  z2([<]),  s(0)  =  c0 

where  [<]  is  the  greatest-integer  function,  has  the  property  that 

x(n  -f  1)  =  4x(n)  —  x2(n),  n  =  0, 1, _ 

If  we  choose  c0  =  4  sin2(7r/9),  then  the  unique  solution  of  this  difference 
equation  is 

x(n)  =  4  sin2  , 

which  has  period  three.  By  the  well-known  result  [8]  which  states  that 
“period  three  implies  chaos,”  the  solution  of  the  above  differential  equa¬ 
tion  exhibits  chaos.  Furthermore,  the  equation  of  Carvalho  and  Cooke 

x'(t )  =  ax(t)(  1  —  #([t])) 

is  analogous  to  the  famous  logistic  differential  equation,  but  t  in  one 
argument  has  been  replaced  by  [t].  As  a  result,  the  equation  has  solu¬ 
tions  that  display  complicated  dynamics  [2].  It  seems  likely  that  other 
simple  nonlinear  EPCA  may  display  other  interesting  behavior. 

The  numerical  approximation  of  differential  equations  can  give  rise 
to  EPCA  in  a  natural  way,  although  it  is  unusual  to  take  this  point  of 
view.  For  example,  the  simple  Euler  scheme  for  a  differential  equation 
x'(t )  =  f(x(t))  has  the  form  xn+i  —  xn  =  hf(xn),  where  xn  =  x(nh) 
and  h  is  the  step  size.  This  is  equivalent  to  the  EPCA 

x'(t)  =  f(x([t/h]h)). 
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Impulsive  differential  equations  and  loaded  equations  of  control  the¬ 
ory  fit  within  the  general  paradigm  of  EPCA.  Another  potential  ap¬ 
plication  of  EPCA  is  the  stabilization  of  hybrid  control  systems  with 
feedback  delay.  By  a  hybrid  system  we  mean  one  with  a  continuous 
plant  and  with  discrete  (sampled)  controller.  Some  of  these  systems 
may  be  described  by  EPCA  [3]. 

EPCA  have  only  been  researched  for  a  few  years.  In  each  of  the 
areas  —  existence,  asymptotic  behavior,  periodic  and  oscillating  solu¬ 
tions,  approximation,  application  to  control  theory,  biomedical  models, 
and  problems  of  mathematical  physics  —  there  appears  to  be  ample 
opportunity  for  extending  the  known  results. 

2.1.  Boundary- Value  Problems.  The  main  emphasis  in  the  project 
was  the  study  of  partial  differential  equations  (PDE)  with  piecewise 
continuous  delay.  The  first  fundamental  paper  [9]  in  this  direction  ap¬ 
peared  in  1991.  It  has  been  shown  in  [9]  that  these  equations  naturally 
arise  in  the  process  of  approximating  PDE  by  using  piecewise  constant 
arguments.  Thus,  if  in  the  equation 

ut  =  a2uxx  —  bu, 

which  describes  heat  flow  in  a  rod  with  both  diffusion  a2uxx  along  the 
rod  and  heat  loss  (or  gain)  across  the  lateral  sides  of  the  rod,  the  lateral 
heat  change  is  measured  at  discrete  times,  then  we  get  an  equation  with 
piecewise  constant  argument  (EPCA) 

ut(x,t)  =  a2uxx(x,t )  —  bu(x,nh), 
t  €  [nh,(n  +  l)h],  n  =  0,1,... 

where  h  >  0  is  some  constant.  This  equation  can  be  written  in  the 
form 


ut(x,  t )  =  a2uxx(x,  t)  -  bu(x,  [ t/h]h ),  (1) 

where  [•]  designates  the  greatest-integer  function. 

The  diffusion-convection  equation 


-  rux 
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describes,  for  instance,  the  concentration  u(x,t )  of  a  pollutant  car¬ 
ried  along  in  a  stream  moving  with  velocity  r.  The  term  a2uxx  is  the 
diffusion  contribution  and  —rux  is  the  convection  component.  If  the 
convection  part  is  measured  at  discrete  times  n/i,  the  process  results  in 
the  equation 


ut(x ,  t. )  =  a2uxx(x,  t )  -  rux(x,  [ t/h]h ).  (2) 

These  examples  indicate  at  the  considerable  potential  of  EPCA  as  an 
analytical  and  computational  tool  in  solving  some  complicated  prob¬ 
lems  of  mathematical  physics.  Therefore,  it  is  important  to  investigate 
boundary- value  problems  (BVP)  and  initial-value  problems  (IVP)  for 
EPCA  in  partial  derivatives,  and  explore  the  influence  of  certain  dis¬ 
continuous  delays  on  the  behavior  of  solutions  to  some  typical  problems 
of  mathematical  physics. 

The  topic  of  [9]  is  the  BVP  consisting  of  the  equation 


du(x ,  t) 
dt 


-f~ 


p  (£) u(i>()=Q  (£Mx’ KW’  (3) 


where  P  and  Q  are  polynomials  of  the  highest  degree  m  with  coefficients 
that  may  depend  only  on  x ,  the  boundary  conditions 


m 

LjU  =  £  1))  =  0,  (4) 

k=l 

(Mjk  and  Njk  are  constants,  j  =  1, . . . ,  m) 
and  the  initial  condition 


m(x,0)  =  uq(x).  (5) 

Here  [•]  designates  the  greatest-integer  function,  ( x,t )  E  [0,1]  x  [0,oo), 
and  h  =  const.  >  0.  Conditions  (4)  will  be  written  briefly  as 

Lu  =  0. 

An  important  result  has  been  established  that  BVP  (3),  (4),  (5)  has  a 
solution  in  [0, 1]  x  [n/t,  ( n  -f  l)/i],  if  the  following  hypotheses  hold  true: 
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(i)  The  boundary- value  problem 


P  (J-)  X  -  XX  =  0,  LX  =  0 

is  self-adjoint,  all  its  eigenvalues  A j  are  positive. 

(ii)  For  each  A  j,  the  roots  of  the  equation  P(z )  —  \j  =  0  have  non¬ 
positive  real  parts. 

(iii)  The  initial  function  uq(x)  £  Cm[0, 1]  satisfies  (4). 

The  solution  un(x,t)  of  BVP  (3),  (4),  (5)  on  the  interval  nh  <  t  < 
(n  -f  1  )h  is  represented  in  the  form  of  a  Fourier  series 


un(x,  t )  =  £  Xj(x)Tnj(t ),  (6) 

j=i 

where  Xj(x)  are  the  eigenvalues  of  the  operator  P.  The  functions  Tnj(t ) 
are  solutions  of  ordinary  EPCA  that  arise  after  separation  of  variables. 

For  instance,  in  [0, 1]  x  [ nh ,  (n  +  1  )h\,  the  solution  un(x,  t )  of  Eq.  (1) 
with  boundary  conditions  un(x,nh)  =  un(x)  is  sought  in  form  (6). 
Separation  of  variables  produces 

Xj(x)  =  y/2  sm(7rjx),  T'nj{t )  +  a27r2j2Tnj(t)  =  -bTnj(nh), 

whence 

Tni(t)  =  Cnje-’‘V^1-^  -  -P^Tnj(nh). 

CL  7T  J 

We  put  t  =  nh  in  this  equation  and  get 

Cnj  =  ^1  +  a2n2j2)  T”j(nh)i 

that  is, 


where 


Tnj(t)  =  Ej(t  -  nh)Tnj(nh), 


b 

Cp"lPj2 


At  t  —  (n  -j-  l)h  we  have 


Tnj((n  +  1  )h)  -  Ej(h)Tnj{nh) 


(7) 
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and  since 

T, ij((n  + 1  )h)  =  T„+ij((n  +  l)ft), 

then 

Tn+ij((n  +  l)h)  =  Ej(h)Tnj(nh) 

and 

Tnj(nh)  =  E?(h)T0j(  0). 

Therefore, 

Tnj(t)  =  Ej(t  -  nh)E?(h)Toj(0) 

and 

un(x,  t)  =  52  V2Ej(h)T0j(0)Ej(t  -  nh )  sin(7rja:).  (8) 

i- 1 

Putting  t  =  0,  n 

=  0  gives 

u0(x)  =  53  Toj(0)v^sin(7r7a;)  dz 
i=i 

and 

7q j(0)  =  V2 u0{x)sin(TTjx)dx. 

If  \Ej(h) |  <  1,  then  solution  (8)  decays  exponentially  as  t  — ►  oo,  uni¬ 
formly  with  respect  to  x.  From  (7)  it  follows  that  this  is  true  if 

pa?ir*h  .  i 

2  2^i^  2  2e  T-  J- 

— a  7T  <  b  <  an — rn — 7- 

ea  *  h  -  1 

Furthermore,  from  the  equations 

Tnj(nh )  = 

=  Ej(h)Toj(0),  T„j((n  +  1)A)  =  E?+'(h)Toj(  0) 

we  see  that  T„j(n/i)Tn,((n-|-l)/i)  <  0  if  Ej(h )  <  0.  The  latter  inequality 

holds  true  if 

a27r2 

^  >  ea***h 

(9) 
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Hence,  under  condition  (9)  each  function  Tnj(t)(j  =  1,2,...)  has  a  zero 
in  the  interval  [nh,  (n  +  l)/i],  in  sharp  contrast  to  the  functions  Tj(t)  in 
the  Fourier  expansion  for  the  solution  of  the  equation  =  a2uxx  —  bu 
without  time  delay.  Moreover,  the  inequality  Ej{h)  <  0  takes  place  for 
sufficiently  large  j  and  any  b  >  0.  Therefore,  for  6  >  0  and  sufficiently 
large  j,  the  functions  Tnj(t )  are  oscillatory. 

Eq.  (2)  on  nh  <  t  <  (n  +  1  )h  becomes 

dun(x,t)  _  2 d2un(x,t) 
dt  a  dx 2 

and  we  differentiate  the  latter  with  respect  to  t  to  obtain  the  equation 

dyn  _  2  d2yn  _  dun 

dt  ~a  dx 2 ’  Vn~  dt' 

whose  solution  is  sought  in  form  (6).  Separation  of  variables  leads  to 
the  equations 

X"(x)  +  XX{x)  =  0,  T'n(t)  +  a2\Tn(t)  =  0, 
and  the  boundary  conditions  un(0,t)  =  un(l,t)  =  0  give  Xj  =  j2 n2  and 


-  ru'Jx), 


yn(x,t)  =  £  V2Tnj(nh)e  °2*2j2(*  "^su^tt jx). 
i=  i 

Since 

yn(x,  nh)  =  a2u"(x)  -  ru'n(x),  un(x)  =  un(x,  nh), 

then 


00  /— 

a2u"(x)  -  ru'n(x)  =  y/2Tnj(nh)  sin(7rjx) 

j= i 

and 

Tnj(nh)  =  -a2-K2j2\/2  Jq  un(x)  sin^jx)  dx 

4-  rnjy/2  Jq  un(x)  cos(7 rjx)  dx. 
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Finally 


un(x,t)  =  un(x)  +  £ 

i= i 


V2Tnj(nh)  (l  -  e-«a»aia(*-"ft))  sinQrjj) 

a2n2j2 


Given  the  initial  function  m(x,0)  =  uo(:r),  we  can  find  the  coefficients 
Toj(O)  and  the  solution  uo(x,t)  on  0  <  t  <  h.  Since  uq(x,Ii)  =  u\(x), 
we  can  calculate  the  coefficients  T\j(h)  and  the  solution  Ui(x,t)  on 
h  <  t  <  2h.  By  the  method  of  steps  the  solution  can  be  extended  to 
any  interval  [ nh ,  (n  +  l)h]. 

The  equation 


iq 


du(x,t) 

dt 


r 

2  mo 


+ >'<*>“  HiW 


is  a  piecewise  constant  analogue  of  the  one-dimensional  Schrodinger 
equation 


If  u(x,t)  satisfies  conditions  (4)  and  (5),  with  m  =  2,  then  separation 
of  variables  produces  a  formal  solution 


«,(*,<)  =  E  C^e-^-^Xjix)  +  P~'Qu„(x), 
j=l 

for  nh  <  t  <  (n  +  1  )h.  Here,  Xj{x)  are  the  eigenfunctions  of  the 
operator  q2(d?/dx2)/2mo ,  and  P~lQun(x)  is  the  solution  vn(x)  of  the 
equation 

q2v^{x)  =  2  moV(x)un(x) 

that  satisfies  (4). 

The  Fourier  method  was  also  used  to  find  weak  solutions  of  the 
boundary- value  problem  (3),  (4),  (5)  and  it  is  easily  generalized  to  simi¬ 
lar  problems  in  Hilbert  space.  First,  we  remind  a  few  well  known  defini¬ 
tions.  Let  H  be  a  Hilbert  space  and  let  P  be  a  linear  operator  in  H  (ad¬ 
ditive  and  homogeneous  but,  possibly,  unbounded)  whose  domain  D(P) 
is  dense  in  H ,  that  is  V(P)  =  H.  The  operator  P  is  called  symmetric  if 
(Pu,  v)  =  (u,  Pu),  for  any  2)(P).  If  P  is  symmetric,  then  (Pu,  v) 
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is  a  symmetric  bilinear  functional  and  (Pu,  it)  is  a  quadratic  form.  A 
symmetric  operator  P  is  called  positive  if  (Pit,  it)  >  0  and  (Pit,  it)  =  0 
if  and  only  if  it  =  0.  A  symmetric  operator  P  is  called  positive  definite 
if  there  exists  a  constant  j2  >  0  such  that  (Pit,  it)  >  72||it||2.  With  ev¬ 
ery  positive  operator  P  a  certain  Hilbert  space  Hp  can  be  associated, 
which  is  called  the  energy  space  of  P.  It  is  the  completion  of  D(P), 
with  the  inner  product  (u,v)p  =  ( Pu,v );  it,  v  G  'D(P).  This  product 
induces  a  new  norm  |jif]|p  =  (Pit, it)1/2,  it  G  T>(P),  aad  if  P  is  positive 
definite,  then  ||tt||  <  7-1||u||p.  Since  D(P)  is  dense  in  H ,  it  follows 
by  using  the  latter  inequality  that  the  energy  space  Hp  of  a  positive 
definite  operator  P  is  dense  in  the  original  space  H. 

Assuming  P  is  positive  definite,  we  may  consider  the  solution  u(x,  t ) 
of  (3),  (4),  (5)  for  a  fixed  t  as  an  element  of  Hp.  If  D(Q)  C  H,  then 
Qu(x,  [t/h]h)  may  be  treated  as  an  abstract  function  Qu([t/h\h)  with 
the  values  in  H.  Therefore,  the  given  BVP  is  reduced  to  the  abstract 
Cauchy  problem 

^  +  Pu  =  Qu  ^  ^  h  j  ,  t  >  0,  u|t=0  =  u0  e  H.  (10) 

If  (10)  has  a  solution,  we  multiply  each  term  by  an  arbitrary  function 
g(t)  G  Hp  in  the  sense  of  inner  product  in  H  and  get  on  the  inter¬ 
val  nh  <t<  (n  +  1  )h  the  equation 

+  (u,g)P  =  ( Qu„,g ),  (11) 

where  un  =  u(nh).  Conversely,  if  u  G  C1((n/i,(n  +  1  )h);T>(P))  for  all 
integers  n  >  0  and  satisfies  (11),  then  it  also  satisfies  Eq.  (10).  Indeed, 
if  u  G  T>(P),  then  (ui9)p  —  {Pui9),  (11)  can  be  written 


nh  <t  <  (n  +  1  )h. 


Since  Hp  is  dense  in  P,  then  u(f)  is  a  solution  of  Eq.  (10). 
Definition.  An  abstract  function  u(t) :  [0,  oo)  —>  H  is  called  a  weak 
solution  of  problem  (10)  if  it  satisfies  the  conditions: 
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(i)  u(t)  is  continuous  for  t  >  0  and  strongly  continuously  diferen- 
tiable  for  t  >  0,  with  the  possible  exception  of  the  points  t  —  nh 
where  one-sided  derivatives  exist. 

(ii)  u(t )  is  continuous  for  t  >  0  as  an  abstract  function  with  the  values 
in  Hp  and  satisfies  Eq.  (11)  on  each  interval  nh  <  t  <  (n  +  l)h, 
for  any  function  g(t ) :  [0,  oo)  — ►  Hp. 

(iii)  u(t)  satisfies  the  initial  condition  (10),  that  is, 

lim|K<)  -  «o|U  =  0- 

A  weak  solution  u(t)  is  also  an  ordinary  solution  if  u(t)  G  D(P),  for 
any  t  >  0,  and  w(x,  t)  — ►  uq(x)  as  t  — ►  0  not  only  in  the  norm  of  H  but 
uniformly  as  well.  It  is  said  that  a  symmetric  operator  P  has  a  discrete 
spectrum  if  it  has  an  infinite  sequence  {Aj}  of  eigenvalues  with  a  single 
limit  point  at  infinity  and  a  sequence  {Aj}  of  eigenfunctions  which  is 
complete  in  H.  Suppose  the  operator  P  in  (11)  is  positive  definite 
and  has  a  discrete  spectrum  and  assume  the  existence  of  a  solution 
u(t)  —  u(x,t)  to  Eq.  (11)  with  the  condition  u(0)  =  uq.  On  the  interval 
nh  <  t  <  (n  +  l)h  this  solution  can  be  expanded  into  series  (6),  where 
Tj(t)  =  (u(t),Xj).  To  find  the  coefficients  Tj(t),  we  put  g(t)  =  Xt 
in  (11)  and  since  Xk  does  not  depend  on  t,  then 

(n,Xt)p  =  (Pu,Xk)  =  (u,  PXt)  =  h{u,Xt)  = 

which  leads  to  the  equation 

Pnj(t)  “b  ^ jPnj(t )  =:  {QumXj) 

By  selecting  a  proper  space  H,  a  weak  solution  corresponding  to  condi¬ 
tions  (4)  can  be  constructed.  A  theorem  has  been  stated  in  [10]  that  if 
P  and  Q  are  linear  operators  in  a  Hilbert  space  and  P  is  positive  defi¬ 
nite  with  a  discrete  spectrum,  then  there  exists  a  unique  weak  solution 
of  problem  (10). 
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2.2.  Initial- Value  Problems.  This  topic  has  been  explored  recently 
by  Wiener  and  Debnath  [10].  Eq.  (3)  with  constant  coefficients  and 
initial  condition  (5)  has  been  considered  in  the  domain 


(x,t)  E  H  =  (—00,00)  x  [0, 00). 

Let  un(x,  t )  be  the  solution  of  the  given  problem  on  nh  <  t  <  (n+l)/i, 
then 


Qt  ^  +  Pun(x ,  t)  =  Qun(x), 

(12) 

where 

un(x)  =  un(x,nh). 

(13) 

Write 

un(x,  t)  =  wn(x,  t)  +  vn(x), 

which  gives 

the  equation 

dt  +  Pwn  +  Pvn(x)  =  Qun(x ), 

and  require 

that 

dwn  _ 

at  +p“'"  =  0’ 

(14) 

Pvn(x)  -  Qun(x). 

(15) 

If  vn(x)  is  a 

solution  of  ODE  (15),  then  at  t  =  nh  we  have 

wn(x,nh )  =  un(x )  -  vn(x), 

(16) 

and  it  remains  to  consider  Eq.  (14)  with  initial  condition  (16).  It  is 
well  known  that  the  solution  E(x,t )  of  the  problem 


—  +  Pw  =  0,  w(x,  0)  =  w0(x ),  (17) 

with  wo(x)  =  6(x),  where  6(x)  is  the  Dirac  delta  functional,  is  called 
its  fundamental  solution.  The  solution  of  IVP  (17)  is  given  by  the 
convolution 


w(x ,  t)  =  E(x,  t)  *  wq(x). 


(18) 
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Hence,  the  solution  of  problem  (14),  (16)  can  be  written  a s 

wn(x,t)  =  E(x,t  —  nh)  *  wn(x,nh),  (19) 

and  the  solution  of  (12),  (13)  is 

u„(x ,  t)  =  E{x ,  t-nh)*  ( un(x )  -  vn(x ))  +  vn(x ),  (20) 

(nh  <t<(n  +  1  )h). 

Continuity  of  the  solution  a,t  t  =  (n  +  l)h  implies 

un(x,  ( n  +  1  )h)  =  un+i(rc,  (n  +  1  )h)  =  un+i(x), 

that  is, 


un+i(x)  =  E(x,  h)  *  (un(x)  -  vn(x))  +  vn(x).  (21) 

Formulas  (20),  (21)  successively  determine  the  solution  of  IVP  (3),  (5) 
on  each  interval  nh  <  t  <  (n  +  1  )h.  Indeed,  from  Pvo(x)  =  Quq(x) 
we  find  ^o(^)  and  substitute  both  uq(x)  and  i;o(z)  in  (20)  and  (21)  to 
obtain  uo(x,t)  and  i/i(a:).  Then  we  use  ui(a:)  in  (15)  to  find  v\(x )  and 
substitute  u\(x )  and  vi(z)  in  (20)  and  (21),  which  yields  it\(x,  t)  and 
uz(x ).  Continuing  this  procedure  leads  to  un(x,t ),  the  solution  of  (3) 
on  [nh,(n  +  1  )h].  The  solution  vn(x )  of  (15)  is  defined  to  within  an 
arbitrary  polynomial  q(x)  of  degree  <  m.  Since  q(x)  is  a  solution  of 
Eq.  (17)  with  the  initial  condition  u;(a:,0)  =  q(x),  then  q(x)  =  E(x,t)* 
q(x),  and  q(x)  cancels  in  the  formulas  (20),  (21).  This  proves  that  if 
Eq.  (17)  with  w;(a:,0)  =  uq(x )  has  a  unique  solution  on  t  £  (0,oo), 
then  there  exists  a  unique  solution  of  IVP  (3),  (5)  on  (0,oo)  and  it  is 
given  by  (20),  for  each  interval  nh  <  t  <  (n  -{-  l)h.  Thus,  there  exist 
unique  solutions  of  Eqs.  (1)  and  (2),  with  u(x,0)  —  uq(x ),  in  the  class 
of  functions  that  grow  to  infinity  slower  than  exp(a:2)  as  |z|  — >•  oo.  For 
Eqs.  (1)  and  (2)  we  have 

vn(x)  =  a~2b  f  (x  -  s)un(s)ds  and  u„(x)  =  a~2r  j  un(s)ds, 

J  0  J 0 

respectively,  and  E(x,t )  =  exp(—x2/4a2t)/2ay/iri. 

Formula  (20)  for  the  solution  of  the  problem 
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on  nh  <  t  <  (n  4-  \)h  becomes 


un(x,  t)  -  ^1  -  E(x,  t-nh)*  un(x)  +  ~^un(x), 


where  E(x,t)  is  the  same  as  in  Eqs.  (1)  and  (2). 

The  above  method  may  also  be  used  to  solve  IVP  for  PDE  of  any 
order  in  t  with  piecewise  constant  delay  or  systems  of  such  equations.  In 
the  latter  case,  P  and  Q  in  (3)  are  square  matrices  of  linear  differential 
operators  and  u(x,t)  is  a  vector  function.  Thus,  the  solution  un(x,t) 
of  the  problem 


uit(x,t)  =  a2uxx(x,t)  -  buxx(x ,  [*]), 
u(x,  0)  =  fo(x),  ut(x ,  0)  =  g0(x) 


on  n  <  f  <  n  +  1  is  sought  in  the  form  u„(x,t )  —  wn(x,t)  4-  vn(x ) 
whence  a2v”(x)  —  6w"(x,n)  =  0  and  d2wnfdt 2  =  a2d2wn/dx2.  Setting 
u(x,n)  =  fn(x),  ut(x,n )  =  gn(x)  gives 

vn(x )  =  a~2bfn(x),  w(x, n)  =  (1  -  a~2b)fn(x),  wt(x, n)  =  gn(x), 


and 

un(x,t) 


-a2fn(x)+(l 


b  \  fn(x  -  a(t  -  n ))  -t-  fn(x  4-  a(t  -  n)) 

)  2 

1  /■*+«(<-«) 


Putting  t  =  n  4-1  produces  the  recursion  relations 

fn+i(x) 


=  +  ('--2 


b  \  fn(x  -  a)  4-  fn{x  +  a) 


) 


g„+i(x)  =  (l  -  A) 


afnjx  +  a)  -  afh(x  -  a) 


2 


1 
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Loaded  partial  differential  equations  have  properties  similar  to  those 
of  equations  with  piecewise  constant  delay.  The  IVP  for  the  following 
class  of  loaded  equations 


du(x,  t ) 
dt 


=p(t)u(x'i)+£Qi{ic)u{x't’)’ 

u(x,  0)  =  u0(x) 


(22) 


was  considered  in  [1]  and  [10],  where  ( x,t )  G  W1  x  [0,T],  the  tj  G 
(0,T]  are  given,  P(s )  and  Q(s)  are  polynomials  in  s  =  (si,...,s„), 
and  E  |Qj(s)|  ^  0.  Eq.  (22)  arises  in  solving  certain  inverse  problems 
for  systems  with  elements  concentrated  at  specific  moments  of  time. 
The  Fourier  transform  U(s,t)  of  u(x,t )  satisfies  the  equation 


Ut(s,t)  =  P(is)U(s,t)  +  £  Qj(is)U(s,tj): 

j=  i 

whence, 


U{s,t)  =  U0(s)ep(“>‘  +  k(P(is),t)  t,  Qj(is)U(s,tj),  (23) 

j=l 

where  Uq(s )  is  the  Fourier  transform  of  uq(x)  and 

k(P(is),t)  =  j‘ep^dy. 

Denote 


Aj  =  U<,(s)epW>,  kj  =  k(P(is),  tj),  B  =  £  Qj(is)U(s,  tj), 

>='  (24) 


then  multiply  by  Qj(is )  each  of  the  equations 

U(s,tj)  —  Aj  +  kjB,  j  =  1, . . . ,  q 
and  add  them.  Hence, 


B  =  i:AjQj(is)  +  BikjQj(is) 
j=  l  i=i 
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or 

(l  -  E  B  =  t  AjQj(i>).  (25) 

V  j= i  /  j=i 

The  equation 

A (s)  =  1  -  £  tj )  =  0  (26) 

;= i 

is  called  the  characteristic  equation  for  (22)  and  its  solution  set  Z  is 
called  the  characteristic  variety  of  (22).  It  is  said  [1]  that  (22)  is  abso¬ 
lutely  nondegenerate  if  Z  =  0,  nondegenerate  of  type  a  if 

a  =  inf  |  Im  s\  <  oo,  s  E  Z  ^  C", 

and  degenerate  if  Z  =  Cn.  The  case  Z  =  0  implies  A (s)  =  const., 
since  A(s)  is  meromorphic,  and  a  meromorphic  function  that  is  not 
constant  assumes  every  complex  value  with  at  most  two  exceptions. 
The  equation  A(s)  =  C  can  be  written  as 

P(is)  +  E  Qj(U)  -  E  Qj(is)epW‘  =  CP(is) 

j= 1  j=l 

and  is  possible  for  q  >  1  only  if  P(s )  =  const,  otherwise  exp (P(is)tj) 
would  grow  faster  than  any  polynomial,  which  breaks  the  latter  equal¬ 
ity.  For  q  =  1  we  have 

.  ,  ,  P(is)  +  Ql(is)  -  Qi(is)ep(i‘)‘‘ 
n{  )  ~  P(is) 

and  in  this  case  Z  =  0  is  equivalent  to  P(is )  +  Qi(^)  =  0.  On  the 
other  hand,  A(s)  =  0  is  equivalent  to 

P(is)  +  E  QiM  -  E  s  0, 

}  =  1  J  =  l 

which  implies  P(s)  =  const.  This  establishes  the  following  proposition 
which  was  stated  in  [1]  without  proof,  namely  Eq.  (22)  is  absolutely 
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nondegenerate  if  and  only  if  either  of  the  following  conditions  holds 
true: 

(i)  P(s)  =  Cu  Qj(s)k(Ci,tj)  =  C2  7^  1; 

3= 1 
or 

(ii)  q  =  1,  P(s)  +  Qx(s)  =  0. 

Eq.  (22)  is  degenerate  if  and  only  if 

P(s)  =  Cu  'tQj(s)k(Cl,ij)  =  l. 
j=  1 

Substituting  B  from  (25)  in  (23)  leads  to  the  proof  that  the  uniqueness 
classes  for  the  solution  of  the  Cauchy  problem  for  an  absolutely  nonde¬ 
generate  equation  (22)  are  the  same  as  those  for  the  equation  (without 
“loads”)  ut(x,t )  =  Pu(x,t).  The  homogeneous  degenerate  IVP  (22) 
(mo(^)  =  0)  has  nontrivial  solutions,  with  compact  support.  Suppose 
that  (22)  is  of  finite  type  a(0  <  a  <  oo)  and  that  u(x,t )  is  a  solution 
of  (22)  with  uq(x )  =  0.  If 

|u(x,<)|  <Ce“w,  ter,  «6[0,T],  (27) 

and  a  <  a,  then  u(x,t )  =  0.  For  any  a  >  a  there  exists  a  solution 
u(x,t )  ^  0  of  (22)  with  uq(x)  =  0  satisfying  (27).  Integral  transforma¬ 
tions  have  also  been  used  in  the  study  of  EPCA. 

Consider  the  nonlinear  initial-value  problem 

^  =  A(D)u{x,t)  +  f{t,u(x,  [*])), 
u(x,0)  =  u0(x), 

where  u(x, t)  and  uq(x)  are  m- vectors,  x  =  (x\,X2 ,xpt)  €  R  , 

A(D )  =  £  A„D° , 

\a\<r 

a  =  (ai,02,...,ajv),  |a|  =  a!+a2H - kaN, 

Da  =  £>?'  •  •  •  D%H ,  Dt  =  id/dxk  (k  =  1,2, . . .  ,N), 


(28) 
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the  coefficients  Aa  axe  given  constant  matrices  of  order  m  x  m,  and 

the  m-vector  /  G  C1([n,n  +  1)  x  £2(RN),£2(RN)),  n  =  0,1,2, _ 

The  number  r  is  called  the  order  of  the  system.  It  is  assumed  that 
wo  E  £2(Rn),  and  the  solutions  sought  are  such  that  u(x,t)  E  £2(RN), 
for  every  t  >  0.  Let  /q (s),/z2(s), . . . , be  the  eigenvalues  of  the 
matrix  A(s).  The  system 

%  =  A(D)n  (29) 

is  said  to  be  parabolic  by  Shilov  if 

Refij(s)  < -c\s\h +  b,  j  =  l,...,m 

where  h  >  0,  c  >  0,  and  6  are  constants.  For  a  fixed  t  we  may  consider 
the  solution  u(x,t )  as  an  element  of  £2(RN),  and  f(t,u(x,[t]))  may 
be  treated  as  an  abstract  function  f(t,u([t]))  with  the  values  in  £ 2. 
Therefore,  IVP  (28)  is  reduced  to  the  abstract  Cauchy  problem 

cfox 

=  Au  +  fit ,  «([*])),  u|<=0  =  M0  G  £2.  (30) 

Applying  to  (29),  with  the  initial  condition  u(x,0)  =  uq(x),  the  Fourier 
transformation  T  in  x  produces  the  system  of  ordinary  differential  equa¬ 
tions 


Ut(s,t)  =  A(s)U(s,t ),  (31) 

with  the  initial  condition  U(s,  0)  =  Uo(s),  where  U(s,t )  =  J^xJ)), 
Uq(s)  =  3r(wo(x)),  and  A(s)  is  a  matrix  with  polynomial  entries  de¬ 
pending  on  s  =  ($i,S2,  •  •  • ,  sn)-  The  solution  of  (31)  is  given  by  the 
formula 


U(s,  t)  =  etA^Uo(s). 

Parabolicity  of  (29)  by  Shilov  implies  that  the  semigroup  T(t )  of  op¬ 
erators  of  multiplication  by  eiA^s\  for  t  >  0,  is  an  infinitely  smooth 
semigroup  of  operators  bounded  in  £2(RN).  Together  with  the  require¬ 
ment  h  =  r,  this  ensures  that  the  Cauchy  problem  for  (29)  is  uniformly 
correct  in  ^(R^)  and  all  its  solutions  are  infinitely  smooth  functions 
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of  tf,  for  t  >  0.  Since  /  is  continuously  differentiable,  problem  (28)  has 
on  [0, 1)  a  unique  solution 

u{t)  =  T(t)u0  +  JQT(t  -  s)/(s,  m0)  ds. 

Denoting  ui  =  u(l),  we  can  find  the  solution 

u(t )  =  T(t  —  1)mi  -f  J |  T(t  —  s)f(s,  Mi)  ds 
of  (28)  on  [1,2)  and  continue  this  procedure  successively.  If 

f{t,u([t]))  =  Bu([t]), 

where  B  is  a  constant  matrix,  the  solution  of  (28)  for  t  £  [0,  oo)  is  given 

by 

<■(*)=  (T(t  -  [t])  +  j‘}T(t  -  s)B ds} 

x  JI  (t(  1)  +  J*_t  T(k  -  s)B  ds)  uo, 

This  proves  that  problem  (28)  has  a  unique  solution  on  IR^  x  [0,oo) 
if  system  (29)  is  parabolic  by  Shilov,  the  index  of  parabolicity  h  co¬ 
incides  with  its  order  r,  and  /  £  +  1])  x  £2(RN),£  (R^)), 

n  =  0, 1, 2, . . . . 

2.3.  Wave  Equations  with  Discontinuous  Time  Delay.  The  m 
fluence  of  terms  with  piecewise  constant  time  on  the  behavior  of  the 
solutions,  especially  their  oscillatory  properties,  of  the  wave  equation 
was  initiated  in  1991  by  Wiener  and  Debnath  [11,  12]. 

First,  we  shall  discuss  separation  of  variables  in  systems  of  PDE. 
Consider  the  BVP  consisting  of  the  equation 

Ut(x,  t )  =  AUxx(x,  t )  +  BUxx(x,  [*]),  (32) 

the  boundary  conditions 

t/(0,  t)  =  C/(l,0  =  0,  (33) 

and  the  initial  condition 


U(x,0)  =  UQ(x). 


(34) 
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Here,  U(x,t )  and  Uq{x)  are  real  m  x  m  matrices,  A  and  B  axe  real 
constant  m  x  m  matrices  and  [•]  denotes  the  greatest-integer  function. 
Looking  for  a  solution  in  the  form 

U(x,t)=T(t)X(x)  (35) 


gives 

T(t)X(x)  =  AT(t)X"(x)  +  BT([t\)X"(x), 

whence 


(AT(t)  +  BT([t]))-lT'(t)  =  X"(x)X~l(x)  =  -P2, 

which  generates  the  BVP 

X"(x)  +  P2X(x)  =  0,  (36) 

X(0)  =  JA(1)  =  0 

and  the  equation  with  piecewise  constant  argument 

T'(t)  =  - AT(t)P 2  -  BT{[t])P2.  (37) 

The  general  solution  of  Eq.  (36)  is 

X(x)  =  cos(xP)C\  -f  sin(xP)C2, 


where 


,  nl  ~(-l)"x2"P2" 
COS (xP)  =  E  - 7^ - . 

n= o  \2n)\ 


.  .  m  ”  (-l)"x2"+1P2"+1 

sin(xP)  =  E - f9„  .  - 

n=o  (2n  +  l)! 


and  Ci,C2  are  arbitrary  constant  matrices.  From  A"(0)  =  0  we  con¬ 
clude  that  C i  =  0,  and  the  condition  X(l)  =  0  enables  us  to  choose 
sin  P  =  0  (although  this  is  not  the  necessary  consequence  of  the  equa¬ 
tion  (sinP)C2  =  0).  This  can  be  written  e‘p  —  e~tP  =  0,  e2tP  =  I.  As¬ 
suming  that  all  eigenvalues  pi,P2, . . .  ,pm  of  P  are  distinct  and  S~l PS  = 
£>  =  diag(pi,p2,  •  •  •  ,Pm),  we  have  exp(2z5D5_1)  =  /,  Se2tT)S~l  =  /, 
and  e2lT>  =  I.  Therefore,  D  —  diag(7r/i,  7172,  •  •  • ,  7r/ra),  where  the  jt  are 
integers,  and  P  =  SlDS-1, 


P2  =  S&S~l 


=  Sdiag(7r2j?,7r% 
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sin(xP)  =  5sin(a:2))5  1  =  S diag(sin 71712;, . . . , sinn jmx)S  l.  Further¬ 
more,  we  can  put 

Pj  =  diag(7r(m(j  -  1)  +  1), ... ,  nmj),  (38) 

0  =  1,2,...) 

in  (36)  and  obtain  the  following  result: 

There  exists  an  infinite  sequence  of  matrix  eigenfunctions  for  B  VP  (36) , 

Xj(x)  =  V2  diag(sin7r(m(jf  —  1)  +  l)x, . . .  ,sin7rmjx),  (39) 

0  =  1,2,...) 

which  is  complete  and  orthonormal  in  the  space  £2[0, 1]  of  m  x  m 
matrices,  that  is, 

J^Xj(x)Xk(x)ds=^  j.*kk 

where  I  is  the  identity  matrix. 

Let  E(t )  be  the  solution  of  the  problem 

T’{t)  =  —AT(t)P2,  T{  0)  =  I  (40) 

and  let 

M(t)  =  E(t )  +  (E(t)  -  I)A~lB.  (41) 

If  the  matrix  A  is  nonsingular,  then  Eq.  (37)  with  the  initial  condition 
T(0)  =  Co  has  on  [0, 00)  a  unique  solution 

T(t)  =  M(t  -  [t])M[t](l)C0.  (42) 

If  ||M(1)||  <  1,  then  ||T’(f)||  exponentially  tends  to  zero  as  t  — >  +00. 
For  the  scalar  parabolic  equation 

ut(x,t.)  =  a2uxx(r,  t )  -t  bu.xx(x,[t]) 

we  have  m  =  1  and  Pj  =  nj,  according  to  (38).  For  Eq.  (40)  with 
A  —  a2  and  P  =  Pj ,  we  have  Ej(t )  =  e  and 
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Hence,  the  inequality  \Mj(l)\  <  1  is  equivalent  to 


■1  <  e-^’  -  4(1  -  e-’*V)  <  1, 

a* 


whence 


—a2  <  b  <  a 


,1 

1  _  e-a***p 


Since  the  function  (1  +  e  *)/(l  —  e  l)  is  decreasing,  all  functions  Tj(t ) 
exponentially  tend  to  zero  as  t  — ►  oo  if  and  only  if 

—  a2  <  b  <  a2.  (43) 

If  b  <  —a2,  then  all  Tj(t)  monotonically  tend  to  infinity  as  t  — >  oo;  and 
if 


b  >  a 


» 1  +  e-®2*2 


1  -  e“a 


2,2  > 


then  all  Tj(t)  are  unbounded  and  oscillatory.  For  any  b  >  a2,  there 
exists  a  positive  integer  jo  such  that  the  Tj(t)  are  unbounded  and  oscil¬ 
latory,  for  j  >  jo.  Indeed,  letting  b  =  a2  +  e  and  solving  the  inequality 


gives 


2  ,  l  +  e-a2*2>2 

a  +€>a  !  _  e-aW 


e-v>'  < 


2 a2  -|-  e’ 

which  holds  for  any  positive  e  and  sufficiently  large  j  and  implies 
Mj(  1)  <  -1.  If  6  =  -a2,  then  M;(f)  =  1,  Tj(t)  =  Tj{ 0),  and 
u(x,t)  =  uq(x),  for  all  t.  Therefore,  the  condition  |6|  <  a 2  is  neces¬ 
sary  and  sufficient  for  the  series 


OO 


u(x,t)  =  £  Tj(t)Xj(x) 
j=  i 


(44) 


to  be  a  solution  of  the  scalar  BVP  (32)-(34),  with  A  =  a2  and  B  =  b,  if 
uq(x)  is  three  times  continuously  differentiable.  The  coefficients  7}( 0) 
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are  given  by 

Tj(0)  =  £  uQ(x)Xj(x)  dx, 

where  Xj(x )  =  y/2sm(irjx)  and  uo(z)  £  C3[0, 1]  satisfies 

w0(0)  =  Uq(  1)  =  0. 

The  solution  T  =  0  of  Eq.  (37)  is  globally  asymptotically  stable  as 
t  — ►  +oo  if  and  only  if  the  eigenvalues  Ar  of  the  matrix  M(  1)  satisfy 
the  inequalities 

|Ar|  <  1,  r  =  l,...,m.  (45) 

If  all  eigenvalues  of  A  have  positive  real  parts  and  Uo(x)  6  C3[0, 1], 
||A~1.B||  <  1,  then  BVP  (32)-(34)  has  a  solution  (44).  This  series  and 
all  its  term-by-term  derivatives  converge  uniformly. 

Separation  of  variables  in  the  equation  with  constant  coefficients 

utt(x,  t)  =  a2uxx(x ,  t)  -  buxx{x ,  [t])  (46) 

and  boundary  conditions  (33)  yields  Xj(x)  =  \/2sin(7r jx)  and  leads  to 
the  EPCA 

T?(t)  +  a2TT2j2Tj(t)  =  bn2j2Tj([t]).  (47) 

For  brevity,  omit  the  subindex  j  and  use  the  substitution  T'(t )  =  V(t), 
which  changes  (47)  to  a  vector  EPCA 

w'(t )  =  Aw(t)  -I-  Bw([t}),  (48) 

where  w  =  co \(T,  V )  and 


Eq.  (48)  on  the  interval  n  <  t  <  n  +  1  becomes 

w'(t )  =  Aw(t )  -1-  Bcn ,  cn  =  w(n) 

with  the  solution 

w(t)  =  M(t  —  n)c„, 
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where 


M(t)  =  eAt  +  (eAt  -  I)A~lB.  (49) 

Therefore,  Eq.  (48)  with  the  initial  condition  w{ 0)  =  Co  has  on  [0,  oo) 
a  unique  solution  given  by  the  right-hand  side  of  (42),  where  M{t )  is 
defined  in  (49). 

For  b  <  0,  the  solution  w  =  0  of  Eq.  (48)  is  unstable.  Indeed, 
computations  show  that 

eAt  =  cos(ut)I  -f  u)~l  sin(a;£)A 


and 


eAt  -  I  = 
where  u  =  anj.  Also 

(eAt  -  I)A~lB  = 

Hence, 


/  ♦  i  _i  • 

cosojt  —  1  u  sine ot 


y  —  usmut  cosut—lj 


b(l  —  cosojt) /a2  0 
^  ( bujsinujt)/a 2  0 , 


coswt  +  6a-2(l  —  cosut)  w_1sina;t 
(ba~2  —  l)u  sin  cut  cos  ut 

and 

det  M(  1)  =  1  — -r  -|-  —t  cos  u ). 

a1  a1 

The  condition  6  <  0  implies  detM(l)  >  1  and  shows  that  at  least  one 
of  the  eigenvalues  A  of  M(  1)  satisfies  |A|  >  1.  Therefore,  ||u/(<)||  — ♦  oo 
as  t  — »  -boo,  for  some  initial  vector  cq  ^  0. 

For  b  >  a2,  the  solution  w  =  0  of  Eq.  (48)  is  unstable.  Calculations 
yield 

det(M(l)  —  A/)  =  A2  —  2  (cos  a;  -I-  sin2  A  +  1  — ~  +  cos  a; 

\  az  2  /  az  az 


M(t)  = 
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and  the  expressions  Ai  =  s  +  d,  X2  =  s  —  d  for  the  eigenvalues  Ai,  A2  of 
M(  1),  where 

.  b  .  2  ^  ,2  (  b  i\-2  b2  .  a  LO 

s  =  cos  u)  +  —  sm  — ,  a  =  I  —  —  1  sin  a;  +  —  sin  — . 

a2  2  Va2  /  a4  2 

The  condition  b  >  a2  shows  that  d2  >  0  and  Ai  >  1.  The  latter 
inequality  implies  ||it;(£)||  — *■  00  as  t  — *  +00,  for  some  initial  vector 
c0  7 ^  0. 

The  solution  ic  =  0  of  Eq.  (48)  is  asymptotically  stable  as  t  — +  +00 
if  and  only  if 

0  <  6  <  a2,  (50) 

and  u>  7^  27rn,  n  =  0, 1, 2, ... .  The  condition  d2  <  0,  which  means  that 
the  eigenvalues  of  M(l)  are  complex,  leads  to 

,w  b 2 

C0S  2  >  (2a2  -  6)2’ 

whence 


6  <  a2  ^1  —  tan2  ^  or  6  <  a2  ^1  —  cot2  ‘j'j  . 

Since  |Ai|  =  { A2I  and  detiVf(l)  =  A1A2,  the  inequality  |Aj|  <  1  is 
equivalent  to  detM(l)  <  1,  that  is,  to  b  >  0.  Therefore,  in  the  case  of 
complex  eigenvalues,  a  criterion  for  asymptotic  stability  is 

i)’“2(i-c°t23)- 

The  inequality  d2  >  0  in  the  case  of  distinct  real  eigenvalues  leads  to 

b  >  max  ^a2  ^1  -  tan2  ,  a2  ^1  -  cot2  ^  j  , 

and  the  inequalities  Ai  <  1,  A2  >  — 1  yield  b  <  a2.  Hence,  in  this  case 
a  criterion  of  asymptotic  stability  is 

2j),a2  (l-cot2  I))  <b<a2. 


max 


a2  (l- 


0  <  6  <  max 
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Finally,  if 

b  =  max  ^z2  ^1  -  tan2  j'j  ,a2  ^1  -  cot2  , 

then  d  =  0  and  Ai  =  A2  =  cos  a;  4-  ba~ 2  sin2  w/2,  whence 

cos  a;  <  Ai  <  cos2u;/2 

and  | Ai|  <  1.  According  to  (45),  this  implies  asymptotic  stability  and 
completes  the  proof  of  criterion  (50). 

If  b  =  a2,  then  Ai  =  1,  A2  =  cos  a;,  and  the  solutions  of  (48)  are 
bounded  but  not  asymptotically  stable.  If  u>  =  2n n,  then  Ai  =  A2  =  1, 
which  leads  to  the  existence  of  unbounded  solutions  for  (48).  If  the 
coefficient  a  is  irrational,  then  (50)  is  a  criterion  of  asymptotic  stability 
of  the  solutions  to  (47)  for  all  j,  since  recalling  that  cj  —  ujj  =  <2717, 
we  note  that  the  equality  anj  =  27m  is  impossible  for  any  irrational 
a.  For  any  rational  a,  there  exist  infinitely  many  integers  j  such  that 
the  corresponding  solutions  wj(t)  of  (48)  are  unbounded.  Furthermore, 
each  component  of  every  solution  of  Eq.  (48)  oscillates  if  and  only  if 
either  of  the  following  conditions  holds  true: 

(i)  b  <  max  |a2  ^1  -  tan2  ,  a2  (l  -  cot2  j  , 

(ii)  max  (a2  (l  -  tan2  j)  ,  a2  (l  -  cot2  )  <  6  < 

and  cosw  <  — 

In  conclusion,  it  is  worth  noting  that  the  asymptotic  properties  of 
Eq.  (47)  depend  on  the  algebraic  nature  of  the  coefficient  a.  For  b  <  0, 
all  solutions  of  Eq.  (47)  are  unstable  and  oscillatory;  for  b  >  a2  all 
solutions  of  Eq.  (47)  are  unstable  and  nonoscillatory.  These  two  cases 
hold  true  for  both  rational  and  irrational  values  of  a.  For 

0  <  b  <  max  ^a2  ^1  —  tan2  ^  ,  a2  ^1  —  cot2  , 
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all  solutions  of  (47)  are  asymptotically  stable  and  oscillatory,  provided 
that  u  ^  27m.  However,  for  any  rational  a,  there  exist  infinitely 
many  integers  j  such  that  u>j  =  27rn,  which  leads  to  the  existence 
of  unbounded  solutions  for  (47).  Furthermore,  since  u  =  u>j  =  anj 
the  inequality  cos  a;  <  —1/2  breaks  down  for  infinitely  many  inte¬ 
gers  j.  Therefore,  under  the  above  hypothesis  (ii),  there  are  infinitely 
many  solutions  of  Eq.  (47)  which  are  asymptotically  stable  and  os¬ 
cillatory,  as  well  as  infinitely  many  solutions  which  are  asymptoti¬ 
cally  stable  and  nonoscillatory  ( u  27m).  Also,  for  u>  ^  27m  and 
a2/ 2sin2(u;/2)  <  b  <  a2,  the  solutions  of  (47)  are  asymptotically  stable 
and  nonoscillatory.  Problems  of  this  nature  deserve  further  investiga¬ 
tion. 
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